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Abstract: 

Charge asymmetries in diffractive electroproduction of two mesons are proportional to 
the interference of Pomeron and Odderon exchange amplitudes. We calculate in the 
framework of QCD and in the Born approximation a forward-backward charge asymmetry 
which turns out to be sizable in a kinematical domain accessible to HERA experiments. 
We predict a distinctive dependence of this asymmetry on the invariant mass of the two 
pions. Testing this prediction is a crucial step in the discovery of the QCD-Odderon. 



1 Unite mixte C7644 du CNRS. 



1. Pomeron and Odderon exchanges are the theoretically dominant contributions 
to hadronic cross sections at high energy. They appear on an equal footing in the QCD 
description of hadronic reactions, and in the lowest order approximation they correspond 
to colour singlet exchanges in the t-channel with two and three gluons, respectively. 

The relevance of the Odderon exchange for hadronic reactions was emphasized long 
ago [ffl. In perturbative QCD the Odderon is described by the Bartels-Kwiecinski- 
Praszalowicz (BKP) equation p]. In spite of many attempts to solve the BKP equation, 
its solutions are still known only partially although much progress has been recently 
achieved || ||, [|, ||. Still this appears to be in a strong contrast with the QCD descrip- 
tion of the Pomeron, where solutions of the leading order Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation are well known M. 

The Odderon remains a mistery also from an experimental point of view. On the 
one hand, recent studies of the elastic pp scattering show that one needs the Odderon 
contribution to understand the data in the dip region . On the other hand, the studies 
of reactions which should select only the odderon exchange didn't show any clear signal of 
its importance. In the case of diffractive r^-meson photoproduction, the QCD prediction 
in the Born approximation for the cross section is rather small M, H. The inclusion of 



evolution following from the BKP equation [[L0| leads to an increase of the predicted 
cross section for this process by one order of magnitude. Recent experimental studies at 
HERA of exclusive 7r° photoproduction |Tl| also indicate a very small cross section for 



this process, in disagreement with theoretical predictions based on the stochastic vacuum 



model In all these meson production processes the scattering amplitude describing 



Odderon exchange enters quadratically in the cross section. 



In |14j] it was suggested to study Odderon effects at the amplitude level by means 
of the asymmetries in open charm production. Since the final state quark-antiquark 
pair has no definite charge parity both Pomeron and Odderon exchanges contribute to 
this process. The Odderon amplitude enters linearly in the asymmetries and therefore 
one can hope that Odderon effects can show up more easily. Moreover, the difficulty 
with the understanding of soft processes in QCD calls for studies of Odderon contribu- 
tions in hard processes, such as electroproduction, where factorization properties allow 
to perturbatively calculate a short- distance part of the scattering amplitude. 

In this paper, we propose to take advantage of a number of interesting features of 
the two pion diffractive electroproduction process to search for the QCD-Odderon at 
the amplitude level. Here again the two pion state doesn't have any definite charge 
parity and both Pomeron and Odderon exchanges contribute. The authors of Ref. [|15[ 



suggest to study the charge asymmetry in soft photoproduction of two pions to select the 
interference of the two amplitudes. Our work shares a number of features with this work. 
The originality of our study of the electroproduction process is to work in a perturbative 
QCD framework which enables us to derive well founded predictions in an accessible 
kinematical domain. 
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Figure 1: Kinematics of the electro-production of two pions 



The aim of the present paper is to study the charge asymmetries in the reaction 

e~(p e ) p(p N ) -»• e~(p' e ) n + (p + ) ?r~(p_) p'(p' N ) (1) 

within perturbative QCD, see Fig. 1. This includes the convolution of perturbatively 
calculable hard subprocess with the two non-perturbative inputs : 2-pion generalized 
distribution amplitude (GDA) and Pomeron-Odderon (P/O) proton impact factors. GDA 



17| are just the light-cone wave functions of the two pion system. Contrary to, say, p- 
meson wave functions, they do not require the selection of a particular charge parity and 
so are ideally suited for studies of the P/O interference. In order to justify the use of 
perturbation theory for this process we consider the electroproduction of this system in 
which the hard scale is supplied by the squared mass — Q 2 of the virtual photon, Q 2 being 
of the order of a few GeV : 
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We intend to study the dominant (for asymptotically large Q 2 ) contribution to the 
asymmetries. This corresponds to the process with longitudinally polarized photon. In 
this case the "longitudinal part" of the two pion wave function is selected, which is 
a straightforward generalization of the longitudinal polarization of vector meson. The 
longitudinal polarization vector is enhanced by a Lorentz boost, together with the cross- 
section of longitudinally polarized pion pair production. This is the twist-2 contribution 
of the subprocess 7*P/0 — > qq — > 7r + 7r~, i.e. in the case of the collinear factorization 
for the qq — > 7r + 7r~ transition. At moderate Q 2 it is customary to take into account 
the transverse momentum k± dependence of the meson wave functions which is a model 
way of accounting for higher twist contributions, so that for growing Q 2 the result is 
undistinguishible from the one of standard collinear factorization Q. The relation between 
the initial and final polarization states is due to the fact that for high-energy diffractive 
processes s-channel helicity conservation is quite well satisfied. Since the cross sections 
a i, and <jt for the longitudinal and transverse photon polarization, respectively, are of 
the same order of magnitude at moderate Q 2 |fL6|| , an experimental separation of ctl is 



highly desirable before confronting our predictions with data. 
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We postpone the consideration of /c^-dependence for future work. 
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Figure 2: Feynman diagrams describing ir + ir electro-production in the Born approximation 

Since the transverse polarization of the pion pair is the only source of the amplitude 
dependence on the azimuthal angle of the pions in their cm. frame, the amplitudes 
and cross sections are independent of this angle in our approximation. As a result, the 
transverse charge asymmetry, resulting from the distribution in this angle and discussed 
in P^| , is zero. Due to that restriction, we only study the forward-backward charge 



asymmetry. 

In the present study we calculate the lowest perturbative order contribution to the 
charge asymmetry, i.e. without taking into account the evolution following from the 
BFKL or the BKP equation. Our results should be therefore treated as an estimate of 
the asymmetries. The above mentioned evolutions can be included into the scattering 



amplitudes in a similar way as in Ref. [1C . 



2. The basic object necessary to calculate the charge asymmetry is the scattering 
amplitude for the process with a longitudinal virtual photon (Fig. 2). 

7l(?) p(Pn) 7r+ (P+) 7T~(P-) p'(pn>) ■ (2) 

We introduce a Sudakov representation with the Sudakov momenta pi,p% obeying the 
equation s = 2pi -p 2 > where s is related to the total energy squared of the virtual photon 
- proton system, Q 2 and the proton target mass M as 

(q + VN? = s-Q 2 + M 2 ws , 

we get for the virtual photon momentum : 

9" = Pf - — P%, (3) 

s 

and for the momentum of the two pion system : 

-Q, 2 _i_ -2 

P21T = V 1 - — m + P2 + V2.L = -P2-, ■ ( 4 ) 

The quark (Zi) and antiquark (/ 2 ) momenta inside the loop before forming two pion 
system are parametrized as : 

K = z P » + m2 + ( l + z P^) 2 p , + {l± + zp27t ±r (5) 

zs 
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= ZVt + ^ + ( ^ ^' pg + H_L + * P*r J" (6) 

where 2/ is the relative transverse momentum of the quarks forming the two pion system 
and z = 1 — z, up to small corrections of the order pf^/s. The collinear approximation 
of the factorization procedure of the description of the two pion formation through the 
generalized distribution amplitude leads to the vector I = in the hard amplitude. 
In a similar way as in ([5]), @ we parametrize the momenta of produced pions as 

P+ = CPl + 7 P2 + (RL + CP^r-lT (7) 

P- = CPl + J P 2 + {-P± + C P2n±r (8) 

where 2p is now the relative transverse momentum of the produced pions, C — P2 P+ is the 
fraction of the longitudinal momentum p 2n carried by the produced 7r + , and ( = 1 — (. 
The variable ( is related to the polar decay angle 9 defined in the rest frame of the pion 
pair by 




/3cos0 = 2C-l, (3= x \l ^ • (9) 

Since the " longitudinal part" of the two pion wave function depends only on the angle 9 
and doesn't depend on the azimuthal decay angle <fr (in the same rest frame of the pair) 
we focus on the calculation of forward-backward asymmetries expressed in terms of 9 (see 
below) . 

The nucleon momentum in the initial state is given by 

M 2 

P^N=P2 + —Pl- (10) 

The squared momentum transfer t = r 2 (r M = p^r — <? M ) can written as 
t -r"--# +t . t ^W 2 + Q 2 

S z 

3. It is well known (see e.g. || and references therein) that for large values of s and 
small momentum transfer t the scattering amplitudes can be represented as convolutions 
over the two-dimensional transverse momenta of the t-channel gluons. 

For the Pomeron exchange, which corresponds in the Born approximation to the 
exchange of two gluons in a colour singlet state, the impact representation has the form: 

M P = -is J (2tt) 2 k\ % ( k i,k 2 )-J P {h,k 2 ) (12) 

where Jj, ^ 7T+7T (£ 1; k 2 ) and Jp~* N '(ki, k 2 ) are the impact factors for transition 

7* — > 7r + 7r~ via Pomeron exchange and of the nucleon in initial state N into the nucleon 

in the final state N'. 
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The corresponding representation for the Odderon exchange, i.e. the exchange of 
three gluons in a colour singlet state, is given by the formula 



_2 



/ d 2 k x d 2 k 2 d 2 k 3 £( 2 )(fci + k 2 + k 3 -p 27T ) r — " ' ~ y \- \ ' 
3! J {^fk\k 2 k\ 



Mo = --^i ^ jg— ( i 3 ) 



where w (hi, k 2 , k 3 ) and Jq~* n ' (k\, k 2 , k 2 ) are the impact factors for the transition 

7* — > 7r + 7T _ via Odderon exchange and of the nucleon in initial state N into the nucleon 
in the final state N'. 



The impact factors are calculated by the standard methods, see e.g. Ref. |L8[ and 
references therein. An important aspect of the present study is the inclusion of an 
appropriate two pion distribution amplitude which we now discuss. 

4. The two-pion generalized distribution amplitude (GDA) [[j/j], [19| contains the 
full strong interactions between the two pions. So far no experimental information exists 
on the two-pion GDA. Watson's theorem imposes that the dynamical phases of the two- 
pion GDA are identical to the phase shifts in elastic tttt scattering as long as m 27T is 
below the inelastic threshold. This relation may be used as an input for a model GDA in 
the m27r-region up to 1 GeV. For higher values of m 2n we assume that the phase is still 
approximately equal to the irn phase shift. The Odderon exchange process involves the 
production of a pion pair in the C-even channel which corresponds to even isospin. In the 
numerical studies we will use a simple ansatz [EDI for the isosinglet distribution amplitude 



$ /=0 (z, (, itl 27T ), in a slightly enlarged m 2n range. We only consider the contributions from 
u- and c?-quarks, i.e. we take n/ = 2. 

A crucial point is the choice of the parametrization of the phases in the GDA's. Let us 
discuss first the isosinglet s— and d— wave phase shifts, 5q and 5 2 . Through interference 
effects, the rapid variation of a phase shift leads to a characteristic m 27r - dependence of 
the asymmetry. We use 



$ /=0 (^, C, m 27r ) = 10z(l - z){2z - 1) R n 

t^ e i5oM \BW fo (ml„)\ +(3 2 e*^) \BW h (m 2 27r )\ P 2 (cos9) 



(14) 



with R n = 0.5 and (3 given by Eq. (^). In our studies we fix the shapes of the phase shifts 
5 and S 2 from a fit to data presented in pl| . Let us stress that the exact description 
of the phase shifts is crucial for our analysis. It is well known that the Breit-Wigner 
parametrization of the amplitude is only a reasonable approximation in the vicinity of the 
resonance peaks. In particular the vicinity of the i^^-threshold leads to some theoretical 
uncertainty of our results, although we do not expect a dramatic effect. 
iBWflm^l is the modulus of the Breit-Wigner amplitudes 

2 

mi 

BW f0 (ml)= ft , m /o = 980MeV, T /o = 50 - 100 MeV (15) 

" L f " L 2tt f 

o mi 

SW ft (mL) = — = — , m f = 1275MeV, r f =186MeV. (16) 

m 2 h -m 2 2lT -im h Y h n n 



5 



To take into account of the uncertainty in the fo width, we will present results for the 
two extreme allowed values. 



The Pomeron exchange process involves the production of a pion pair in the C-odd 
state. Its amplitude can be fully computed for values of m 27r where the timelike electro- 
magnetic pion form factor F^m^) is known provided the lowest Gegenbauer polynomial 
component is dominant. The modulus of F n has been well measured in the process 
e + e" — > 7r + 7r~. By Watson's theorem its phase is equal to the p-wave phase shift Si, pro- 
vided that m 27r is in the range where tttt scattering is elastic. This is rather well satisfied 
for values of m 2?r up to 1 GeV. We assume that for m 27r up to 1.5 GeV we can still use 
the phase shifts from tttt elastic scattering as the phase of the distribution amplitude. In 
our numerical studies we take a i^-parameterization inspired by the N = 1 model of Ref. 



22] 



$ /=1 (z,C,m 27r ) = 6z(l - z)PcoseF„{m% n ) , (17) 



where 



with 



2 1 1 + 1.85 ■ IP' 3 ■ BW„ 

< m 2J = M 7^J^ BW P 



(1 - 0.145) 1 1 + 1.85- 10" 3 



m 2 



BW p (ml) = — , ' - , (19) 



, . _ m 2 p K-4m^) 3 / 2 



r,K =r p ^ f , ?' /2 , m p = 773MeV, r p = 145MeV, 
m 27r W-4m;)' ) '' ! 



and 



2 

m 



BW^ml) = — — , m w = 782 MeV , T w = 8.5 MeV . (20) 

The parametrization of the pion form-factor given by Eq. flTSp leads to a reasonable 
description of the data on the square of the pion form-factor, see |22[| . It describes also 
satisfactorily the p-wave phase shift in the region of 2n invariant mass smaller than 1.5 
GeV. Above 1 GeV the phase shift given by the original parametrization of Ref. [p^j 
strongly overestimates the data points presented in [21], due to a problematic p'(1370) — 
contribution. Since a correct treatment of phases is crucial for our predictions, we do not 
include any p'— contribution. 



5. After choosing the two pion distribution amplitude, the calculation of the 
necessary impact factors is straightforward. Skipping unessential details, let us now 
present the final results. For the longitudinal polarization of virtual photon we obtain 
for the impact-factor for — > 7r + 7r~ the following formulas: 

• Pomeron exchange: 



l 

J?^ + *~ = - * %** Q J dzzzPpihM&^Cml) (21) 

o 
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Pp{h,k 2 ) 



Z^ + fi 2 Z^+fl 2 _ z ^)2 + ^2 fix - Zp 27: ) 2 + /X 2 



',22) 



where &i + fc 2 = and /i 2 = m 2 , + z z Q 2 . 
Odderon exchange (where ki + k 2 + k 3 = p 27T ) 



n (ki,k 2 ,h) 



i e g 3 d abc Q 



AN, 



c 



dz zz P (kiMM) -$ J -°(z,C,mU 



(23) 



Po(ki,k 2 ,k 3 ) 



1 



1 



1 1 



i V (fci - zp 2w ) 2 + /j, 2 (h - zp 2w ) 2 + [i 2 



(24) 



The value of the strong coupling constant g in the hard block is assumed to correspond 
to the 1-loop running coupling constant with rif = 2, a s (Q 2 ) = j- = 127r/[291n(^ — )]. 

In order to estimate the theoretical error we vary the value of Aqcd = 0.2 — 0.35 GeV. 
This does not lead to any dramatic change in the magnitude of the asymmetry. 

Finally we have to fix the soft parts of our amplitudes, i.e. the proton impact factors. 
They cannot be calculated within perturbation theory. In our estimates we will use 



phenomenological, eikonal models of these impact factors proposed in Refs. |23| and |24 
We take 



for the Pomeron exchange : 



J 



N-+N' 



% — 7— 3 



2N, 



c 



A 2 



A 2 + 1/2 pi A 2 + l/2(k 2 + k 2 



(25) 



for the Odderon exchange : 



p.3 Jabc 



AN, 



c 



F{p2n, 0,0) -J] F{k u p 2n - ki, 0) + 2 F(k h k 2 , k 3 



(26) 



i=l 



where 



F(ki, k 2 , k 3 



A 2 



A 2 + \ 



(h - k 2 y + (k 2 - k 3 ) 2 + {k 3 - k,y 



(27) 



and A 
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In these formula, we have denoted the QCD-coupling constant as g. While it is 
natural to take Q 2 as the scale of the strong coupling constant in the upper impact 
factor, a typical hadronic scale M 2 is better suited for the lower one. The value of the 
coupling constant g = g{M 2 ) is one of the main sources of theoretical uncertainties of 
our numerical results. 

In the original Refs. [23|, [24j it was assumed that a so f t — f- « 1. In view of the results 
of the recent studies Jl2|], it seems that this value is too large for a correct description 
of the pp differential cross-section in the region of the dip, and one should rather take 
c^soft = 0.3 — 0.7. In order to visualize this rather large uncertainty we present our results 
with an error band corresponding to this interval of a so f t . We also want to emphasize 
that an increase of the value of a so f t to a so f t ~ 1 raises our predictions for the asymmetry 
by a factor « 1.3. 

6. Let us now present our estimates of the charge asymmetry. We define the 
forward - backward asymmetry as 



2 2 , _ f cos9da(s,Q 2 ,t,ml n ,0) 
W " m2j " fda(s,Q 2 ,t,ml,6) 

J cos 6d cos 6 2 Re Mp L (Mo L ) 



which may be rewritten as 



/ dcosO \Mf\ 2 + \M%\ 2 



(28) 



A(Q 2 



(l+/3)/2 

/ dCj(2C-l)2He 

(l-/3)/2 



Ml L {M lL 



o 



1-MoT 



(29) 



(l+/9)/2 
(l-/3)/2 

Instead of weighted integration of the cross-section it is possible to perform a full 
angular analysis. The numerator of the asymmetry would then be provided by the cos 9- 
term which is characteristic of the longitudinal polarization of the pion pair, while in 
the denominator one may extract experimentally the contribution of the longitudinally 
polarized pion pair in complete analogy to the case of longitudinally polarized vector 
mesons. 

We checked that the squared Odderon contribution in the denominator can be ne- 
glected, so that the asymmetry is practically a measure of the ratio of the Odderon and 
the Pomeron amplitudes. 

There is no s— dependence in our framework, within the approximation which we 
make, provided s is large enough for the usual high energy approximation to hold. The 
charge asymmetry is plotted in Figs. 3 and 4 as a function of the two pion invariant 
mass m 27r . The main characteristic is the high numerical value that we get for values of 
mi-K around the fo mass. This is in strong contrast to the result obtained in a completely 
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Figure 3: Charge asymmetry given by Eq. (|28|) for a minimal /o width, with an error band 
showing the uncertainty comming from different values of a so ft and Aqcd- 



different framework by the authors of Ref. [ 15| . The characterisic m 27r dependence in 
Fig. 3 and Fig. 4 is completely understood in terms of the tttc phase shifts and the factor 
sin(5°— 5 1 ). The phase difference vanishes for m^-K ~ 0.75 GeV and m,2n ~ 1 GeV resulting 
in two zeros of the charge asymmetry. The magnitude of the asymmetry depends much 
on the width of the fo meson which is estimated to be 50 — 100 MeV. The curve shown 
in the Fig. 3 and in the Fig. 4 are obtained with the extreme values of this width. The 
d— wave contribution shows up around 1.3 GeV on Figs. 3, 4 and 5. In Fig. 5 we show 
also the Q 2 dependence of the asymmetry which turns out to be rather moderate. 

The t dependence of the asymmetry, which is plotted in Fig. 6, is quite interesting 
since it has a characteristic zero around t = —0.1 GeV 2 . This zero in the odderon 
amplitude has already been discussed in Ref. flO |. The practical outcome of this t- 



dependence is that it is better to focus on moderate values of t and, if an integration 
over t is needed in order to improve the statistics, to choose at least —t ~ 0.1 GeV 2 as a 
lower integration limit. 

7. Let us now comment on the possible theoretical uncertainties of our calculations. 
• Higher twist contributions and corrections to the calculated expressions of order 



m 



2k 



/Q 2 or t/Q 2 may well be non negligible. For instance, k± effects in the 



7T 1 7T 



wave function mentioned above may be important at relatively low Q 2 |25] as well 
as the effects of transversely polarized photon. An estimate of these corrections is 
certainly desirable but clearly out of the scope of the present paper. 

QCD evolution a la BFKL and \n(s/ Aq CD ) corrections may be calculated for the 
Pomeron and Odderon exchanges as well as the ln(Q 2 /Ag CD ) corrections for the 
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Figure 4: Same as Fig. 3 but for maximal /o width. 

generalized distribution amplitudes. We expect that they do not have drastic effects 
on ratios such as the charge asymmetry which we have calculated. This may be 
controlled in the future. 

The exact values of the scales of the coupling constants should be determined by the 
presently unknown higher order corrections in the hard and soft parts, respectively. 
It is interesting that in the case of the dominance of Abelian contributions to 
radiative corrections, the coupling renormalization would be reduced, by use of the 
Ward identities, to bubble insertions into the gluon line, so that the scales would 
be common for the hard and soft parts and determined by the gluon off-shellness. 
Due to the integration over the gluon momenta the latter should be averaged to 
some intermediate value between Q 2 and M 2 , so that the above naive estimate may 
still be valid. 

The 7T7T distribution amplitude is a non perturbative object which we certainly do 
not exactly know. Although its phase is theoretically under control in the lower 
mass range, its magnitude and its z— and (— dependence may be quite different 
from the simple ansatz that we have adopted. Let us stress however that parts of 
this uncertainty can be resolved by other experiments, namely two pion production 
in e7 collisions [20] and in ep collisions at medium energies [^6 sensitive to the 



C-even and C-odd components of GDA. 

Let us now briefly indicate possible future studies closely related to the approach that 
we have developed. 



Single spin asymmetries should show up in the same reaction ||27|| . Longitudinally 
polarized electron beams are providing circularly polarized virtual photons. In turn, 
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Figure 5: m^-dependence of the asymmetry for t = —.8 GeV 2 and for different values of Q 2 : 
1 GeV 2 (dashed line), 3 GeV 2 (dotted line), 10 GeV 2 (solid line); the d —wave contribution at 
Q 2 = 3 GeV 2 is shown with a dense dotted line; the /o width has been taken to be 75 MeV. 



they give rise to single spin azimuthal asymmetries, similar to those appearing in the 



7*7 production of pion pairs |20| and triplets [^8| . The expression for the asymmetry 
is similar to the one for the charge asymmetry with the notable difference that the 
imaginary, instead of real part of the relevant product of amplitudes appears in 
the numerator of the analog of Eq.^8[ The consequence is that, in the case of the 
Pomeron-Odderon interference, whose amplitudes are mostly imaginary and real, 
respectively, there appears a factor cos(5o — Si) instead of sin(5 — Si), which leads 
to a completely different m 27r dependence. Note that in the case of the calculations 
of the single spin asymmetry one has to take into account transverse polarizations 
of a photon. 

Charge asymmetries can also be studied in the interesting case of 7*7* scattering 
which may be measured at e + e~ colliders. A specific feature of this case is that 
both impact factors are calculable within perturbative QCD. 



8. In conclusion, let us stress that we have demonstrated that the understanding of 
diffractive processes within perturbative QCD is bound to the discovery of sizable charge 
asymmetries in electroproduction of two charged mesons. 

We applied the powerful tool of QCD factorization when the hard part (subprocess) 
is calculated perturbatively, while the soft ingredients (GDA and proton impact factor) 
should be modeled or, better, measured, which poses a new challenging problem for 
exp er iment alist s , 

Data on this reaction in the kinematical domain suitable for our calculation (i.e. large 
s, small t, Q 2 above 1 GeV 2 and m 27r below 1.5 GeV) should be easy to get and analyze 
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Figure 6: t-dependence of the asymmetry for Q 2 = 3 GeV 2 and = 0.95 GeV; the /o width 
has been taken to be 75 MeV. 



by the experimental set-ups HI |TT| and ZEUS |29] at HERA. We are eagerly waiting for 
this confrontation of theory with data, which should lead us towards the discovery of the 
Odderon. 
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